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I.  INTRODUCTION 


Modeling  of  Free  Form  Surfaces 

The  modeling  of  complex  three-dimensional  surfaces  used  in  the 
manufacture  of  cars,  boats,  airplanes  and  other  designs  requiring 
freeform  surfaces  has  traditionally  been  done  using  solid  model  mock- 
ups.  Mock-ups  or  "hard  models"  are  very  costly  to  build  and  difficult  or 
impossible  to  iteratively  modify. 

"Soft  models"  are  computer  data  sets  that  can  be  used  to  mathematically 
generate  freeform  curves  and  surfaces.   In  general  there  are  two  types 
of  representations  used  for  soft  models:  explicit  functional  forms  and 
implicit  functional  forms.   Explicit  functional  forms  are  based  on 
combinations  of  explicit  function  in  which  the  curves  and  surfaces  are 
generated  by  directly  varying  the  geometric  parameters.   Implicit 
functional  forms  are  based  on  parametric  representations  of  the 
explicit  variables.  Parametric  variables  are  varied  to  generate  the 
curves  and  surfaces. 

A  parametric  representation  consists  of  a  set  of  vectors  or  control 
points  and  a  functional  form,  usually  a  polynomial,  in  the  parameter. 
The  control  points  serve  as  weights  to  the  terms  in  the  functional  form. 
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The  first  parametric  representations  used  to  represent  curves  and 
surfaces  were  Bezier  curves  and  surfaces.   This  representation   is  based 
on  a  Bernstein  polynomial  basis  and  is  a  global  approximation  scheme. 
The  number  of  control  points  is  the  same  as  the  order  of  the  curve  and 
changing  any  one  of  the  control  points  affects  the  shape  of  the  entire 
curve. 

A  more  general  extension  of  the  Bernstein  basis  is  the  B-spline  basis,  a 
local  approximation  scheme  which  can  have  an  unlimited  number  of 
control  points.   Only  a  limited  number  of  control  points,  equal  to  the 
order  of  the  curve,  have  an  effect  on  a  given  point.   If  the  number  of 
control  points  is  equal  to  the  order  of  the  curve,  the  B-spline  curve 
reduces  to  the  Bezier  form.   There  are  two  types  of  B-spline  curves  and 
surfaces  that  are  commonly  used  in  computer  aided  design  (CAD),  integral 
and  rational.  The  integral  B-spline  is  a  non-homogenous  vector  valued 
function  whose  components  are  parametric  polynomials.   The  rational 
B-spline  is  a  homogenous  vector  valued  function  whose  components  are 
parametric  polynomials. 

The  first  part  of  this  work  will  deal  with  the  mathematical  basis  for 
the  B-spline  function.   Next,  the  application  of  the  B-spline 
approximation  scheme  to  the  representation  of  geometric  models  will  be 
shown.   Finally,  a  discussion  of  the  program  used  to  generate,  modify, 
and  display  B-spline  surfaces  will  be  presented. 
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Implementations  of  the  Present  Work 

The  surface  modeling  features  that  have  been  implemented  in  this  present 
work  are:  algorithms  to  generate  B-spline  surface  streamlines  and 
normals,  a  primitive  display  file,  an  interactive  environment  to 
globally  change  the  data  set  for  viewing,  and  an  interactive  environment 
to  locally  modify  the  data  set  to  change  the  shape  of  the  surface. 

The  B-spline  algorithms  implemented  are  those  described  by  deBoor  (ref 
1)  for  the  evaluation  of  the  B-spline  approximation  and  its  derivatives. 
The  extension  of  the  algorithms  to  generate  curves  and  surfaces  are 
similar  to  those  described  by  Riesenfeld  (ref  2).   The  display  file 
includes  the  B-spline  control  mesh,  streamlines  in  both  surface 
directions,  and  the  normals  at  the  intersection  of  each  of  the 
streamlines.   The  environment  to  globally  change  the  data  set  for 
viewing  permits  the  user:  1)  to  display  the  data  in  an  isometric,  2) 
rotate  the  data  about  the  Cartesian  axes,  or  to  scale  the  data.   The 
environment  to  locally  change  the  data  is  based  on  two  ideas:  first,  the 
local  support  properties  of  the  B-spline  basis,  and  second,  the  use  of  a 
surface  cursor  rather  than  a  screen  cursor  to  indicate  position  in 
space. 

In  the  next  section  the  mathematical  basis  for  the  B-spline  basis 
function  and  the  form  of  the  B-spline  approximation  scheme  will  be 
presented. 
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II.  B-SPLINE  APPROXIMATION:  DEFINITIONS,  FUNCTION,  PROPERTIES,  AND 

EVALUATIONS 


Curves  and  surfaces  are  generally  expressed  as  functions  of  the  x,  y,  z, 
Cartesian  coordinates.   These  functions  can  be  either  explicit: 

f(x,y,z)  =  0.0 
or  parametric 

f(x)=x(t)  ;  f(y)=y(t)  ;  f(z)  =  z(t). 
These  functions  are  characterized  by  the  degree  or  order  (degree+1)  of 
the  variables.   The  set  of  all  functions  of  the  form 

P  =  a~  +  a.  t.  +  a0t0  +  ...  +  a  t 
n    0    1  1    2  2  n 

is  called  the  polynomials.   The  polynomials  form  a  linear  space  in  the 

mathematical  sense.   They  are  closed  under  the  operations  of  addition 

and  scalar  multiplication.  A  basis  for  the  polynomial  space  will  be  a 

set  of  functions  whose  linear  combination  can  be  used  to  represent  all 

the  members  of  that  space.   There  are  many  different  bases  that  can  be 

used  to  represent  the  polynomial  space.   One  class  of  functions  that  can 

be  used  as  a  basis  for  the  polynomials  is  called  splines.  For  a  set  of 

functions  to  be  a  spline,  they  must  satisfy  the  following  definition 

from  Riesenfeld  (ref  2). 

Definition:  Let  T=  (tQ,t.  ,.  • .  ,t.  )  be  a  vector  of  reals  such  that 

t.<  t.  ..  A  function  N  is  called  a  (polynomial)  spline  function  of 

degree  K-l  (order  K)  if  it  satisfies  the  following  two  conditions: 
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1.  N(K,T)  is  a  polynomial  of  degree  K-l  on  each  sub interval 

(t.,t.+1). 

2.  N(K,T)  and  its  derivatives  of  order  l,2,...,K-2  are 

•   „  •   „(K"2) 
everywhere  continuous,  that  is  N  is  C 

The  points  t.  are  called  knots  and  T  is  called  the  knot  vector.   The 
1 

B-spline  functions  form  a  basis  for  the  polynomials,  and  the  B-spline 

approximation  uses  this  basis  as  its  functional  form.   The  B-spline 

approximation  will  be  used  in  the  geometric  applications  to  represent 

curves  and  surfaces. 


Defining  the  B-spline  approximation's  functional  form. 


The  B-spline  approximation  has  the  following  form: 


F(t)=sumover(i:0,n)  A.  *  N.  „(t) 

l    i,K 

1.  F(t)  :  is  the  parametric  function 

2.  A.  :  is  an  ordered  set  of  constants 

i 

3.  N.  _,(t)  :  are  the  B-spline  basis  functions  N(K,T) 

1  ,K. 

4.  K  :  is  the  order  of  the  function  (degree+1) 

5.  T  :  (t_,t.,  .  .  .  ,t  )  is  the  knot  vector 

6.  n  :  is  the  number  of  constants  (zero  indexed) 
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Details  of  the  B— spline  approximations 

The  use  of  the  B-spline  approximation  will  require  the  evaluation  of  the 
function  and  its  derivatives  with  respect  to  the  parameter.   To  evaluate 
the  B-spline  approximation  requires:  1)  the  order  of  the  basis,  2)  the 
basis  functions,  3)  the  knot  vector,  and  4)  the  constants  A.. 

The  order  of  curve  is  equal  to  the  degree+1  or  the  number  of 
coefficients  in  the  associated  K-l  degree  polynomial. 

The  B-spline  basis  functions  are  defined  by  the  following  recursion 

formula  (ref  1) 

N.  _(t)  =  t  -  t.      N.  _  .(t)  +  t.    -  t    N.  -  _  .(t) 
l,K  _i_       i,K-l       _i*_.        l+l  , K-l 

'i+K-1  "  ti  'i+K  "  ti+l 

N.  ,(t)  -  +  1.0  if  t.<=  t  <  t.  .,  . 


i.  .u;  -  t  1.0  if  t.<= 
1,1      +     elseVo 


i+1 


The  evaluation  of  the  basis  functions  through  order  4  is  given  in 
appendix  A  and  plots  of  the  conical  basis  functions  through  order  4  are 
shown  in  Figure  2-1.   The  basis  functions  are  non-negative  and  the  sum 
of  all  N.  -.(t)  for  any  value  of  t  is  equal  to  1.0   (ref.  2).   These 

1,K. 

properties  are  responsible  for  convexity  of  the  curves  and  surfaces  in 
geometric  applications.   The  definition  given  leads  to  a  one-sided  basis 
function  evaluation,  but  this  does  not  effect  the  continuity  at  the 
point  t=t .  . .   The  number  of  basis  function  for  a  periodic  B-spline  is 

{N.  _}  (i:0,m) 
l  ,fl» 

where  m  is  the  range  of  the  knot  vector  without  multiplicity. 
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The  number  of  basis    functions   for  a  non-periodic   B-spline   is 
{N.   _}    (i:0,K+m-l) 

1,K. 

where  m  is  the  same  as  for  the  periodic  case. 


The  knot  vector  is  a  sequence  of  real  numbers  such  that  t.<=  t.  ,.   If 

1    l+l 

they  are  integers,  the  knot  vector  is  uniform,  otherwise  the  knot  vector 

is  considered  non-uniform»   A  uniform  knot  vector  can  be  used  without 

any  loss  of  generality  in  the  definition  of  the  function  (ref  2).   The 

B-spline  basis  can  be  periodic  or  non-periodic.   If  it  is  periodic  the 

knot  vector  has  the  form 

T:(0,l,2,...m+1) 

and  the  basis  functions  are  shown  in  Figure  2-2.   If  the  knot  vector   is 

non-periodic  then  the  form  is 

T:  (QlQ Q,l,2,...tm-l,Wtffiax t^) 

K  K 

The  basis  functions  are  shown  in  figure  2-3.   The  knot  vector  may  have 

multiplicity  up  to  K,  that  is  t.=t.   =...=t.  _  ..   The  effect  of 

l  l+l      i+K-1 

multiplicity  is  to  reduce  the  differentiability  of  the  basis  functions 

by  one  degree  for  each  multiple  of  t.  as  can  be  seen  in  Figure  2-4,  or 

similarly,  it  reduces  the  order  of  the  curve  on  the  span. 

The  constants  can  be  scalar  valued  or  r  dimensional  vectors.   In 
graphics  applications  the  constants  are  usually  coordinate  vectors  with 
the  function  F(t)  being  the  point  on  a  parametric  curve  or  surface.   The 
characteristic  that  permits  the  efficient  evaluation  of  the  B-spline  is 
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the  local  nature  of  the  basis  function.  That  is,  a  maximum  of  K  basis 

functions  will  be  non-zero  for  any  value  of  t.   In  computing  F(t),  the 

summation  need  only  extend  over  the  K  non-zero  basis  functions  N.  _(t) 

1  ,K. 

(i:j-(K-l),j)  or 

(*) 
F(t)=sumover(i: j-(K-l),j)  A.  *  N.  __(t)    j  s.t.  t .  <=  t  <  t .  .  . 

11  ,K      J       j         j+1 


The  evaluation  of  the  parametric  derivatives  require  a  reduced  set  of 
constants  and  the  lower  order  basis  functions  of  the  span.  The  first 
derivative  is  evaluated  in  the  following  manner 

F(1)(t)  =  (K-l)  sumover(i:0,n))  A^X)  *  N.  _  .(t) 

1      i,K-l 

with  A^1}  =  (A.  -  A.  .)  /  (t.rf  .  -  t.). 
l       l     l-l       i+K-1     l 

In  general  the  derivatives  are  evaluated  as  follows 

F(d^(t)  -  (K-l)*(K-2)*...*(K-d)  sumover(i:0,n)  A^d)  *  N.  ,_  jN(t) 

l      i,(K-d) 

and  the  constants  are 
A<0)  -  A. 

1  1 

.(d)  ,  ,A(d-l)    A(d-1).  .    .  „  ,   .  s  . 

A.   ■  (A.     -  A     )  /  (t>Jjr  ,-  t.)  :d  >  0 
l       l        l-l        l+K-d    l 

with  only  i  =  j-(K-d)+l . . .  j  of  N.  ._  .(t)  being  non-zero  and 

(*)  For  a  particular  t,  j  is  selected  such  that  (s.t.) 
t.  <=  t  <  t.  .  and  the  range  of  the  summation  is  j-(K-l)  (K  terms). 
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F(d)(t)  =  (K-l)*...*(K-d)*  sumover(i:j-(K-d)+l,j)  A^d)  * 


Ni,(K-d)(t) 


with  j  s.t.  t.  <=  t  <  t.  .. 
J         J+l 


The  evaluation  of  the  derivative  exhibits  the  same  properties  as  the 
evaluation  of  the  basis  functions. 


Properties  of  B-spline  functions 

B-spline  functions  exhibit  some  unique  properties  that  are  useful  in 
geometric  applications  such  as  1)  local  support,  2)  (K-2)  continuity  and 
3)  they  are  variation  diminishing.   The  property  of  local  support  is  the 
result  of  the  fact  that  only  K  basis  functions  are  non-zero  for  any 
value  of  the  parameter.   The  (K-2)  continuity  is  guaranteed  from  the 
definition  of  the  general  spline  function.   The  variation  diminishing 
property  (ref  2)  is  described  as  the  ability  of  the  approximation  to 
reproduce  linear  functions  exactly  and  to  have  no  more  zeros  than  the 
primitive  itself. 


Evaluations  of  the  B-spline  Functions 

Algorithms  for  the  efficient  evaluation  of  the  B-spline  and  its 
derivatives  based  on  the  basis  function  definition  were  given  by  deBoor 
(ref  1).   The  algorithm  evaluates  all  the  N.  _(t)  that  are  not  zero  for 
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a  given  t  such  that  t.  <■  t  <  t.  ,.   These  are  shown  in  the  table  2.1 

°  1         l+l 

The  terms  in  the  (K-d)   column  of  this  table  are  the  basis  functions 
needed  to  evaluate  the  d   derivative  of  the  function  F(t)  at  t. 


Ni,l(t)   Ni-l,2(t)   Ni-K+2,K-lCt>   "Mtrtjt(t) 

Ni>2(t)     Bi-K+3,K-l(t)   Ni-K+2,K(t) 


Ni,K-l(t)       Ni-l,R(t) 


Table  2-1.  Required  basis  function  evaluations 

To  generate  the  entries  in  this  table  column  by  column  you  must  first 
define  the  following 

N(r,s)=N.      (t) 
i+r-8 ,  s 

DP(r)=  t.    -  t    + 
l+r        | 

DM(r)  =  t  -  t.  .    |r=l R 

l+l -r  + 

and  the  needed  table  entries  are  N(r,s)  r  ■  l,...,s   s  ■  1  , . . .  ,K 

with  N(r,s)  ■  0.  for  r  >  s  or  r  <  1 . 


-10- 


Substitute  these  into  the  recursion  formula 

N(r,s+1)  -  DM(s+l-r+l)  *  N(r-l,s) 

DPrr-lJ'+'DMU+l-r^i) 

+  DP(r)  N(r,s) 

DP7r)~+~DM(s+f-r7 

The  algorithm  to  generate  the  N(r,s)  entries  in  table  2-1  is  given  in 

Table  2-2  from  deBoor  (ref  1). 

Set  N(l,l)  =  1.0 

for  s  =  1 ,  K-l  do: 

set  DP(s)=t.    -  t,  DM(s)  -  t  -  t .  . 

l+s  i+1-s 

set  N(l,s+1)=0. 

for  r  =  1 , . . .  ,s  do: 

set  M=N(r,s)/(DP(r)  +  DM(s+l-r)) 
set  N(r,s+1)  -  N(r,s+1)  +  DP(r)  *  M 
set  N(r+l,s+l)  =  DM(s+l-r)  *  M 

Table  2-2.  Algorithm  to  evaluate  required  basis  functions. 

Having  calculated  all  the  required  basis  functions,  the  evaluation  of 
the  B-spline  function  and  its  derivatives  requires  the  summation  of  the 
appropriate  products  of  derivative  constants  and  basis  function. 

The  one-sided  definition  of  the  basis  function  generates  an  exception  in 

the  evaluation  of  the  non-periodic  function  when  t  =  t    .At  that 

max 

point   all  the  N.  _(t)  ■  0.  except  N  v(t)   which  is  equal  to  1.0,  and 
the  function  is  equal  to  the  n   constant  at  the  point. 


-11- 


In  the  next  section,  the  B-spline  approximation  will  be  applied  to  the 
modeling  of  freeform  curves  and  surfaces. 
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a.  order  equal  1 


b.  order  equal  2 


c.  order  equal  3.  d.  order  equal  4, 

Figure  2-1  Conical  basis  functions  order  1  thru  4. 
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Figure   2-2        Periodic   basis    functions   of  order  3 
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Figure   2-3        Non-periodic   basis    functions   of  order  3 
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Figure   2-4        Non-periodic  basis    functions   of  order  3  with 
multiplicity. 
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III.    GEOMETRIC  APPLICATIONS  OF  B-SPLINE  APPROXIMATION  TO  COMPUTER 
GENERATED  CURVES  AND   SURFACES 


The  geometric   application     of   the  B-spline   approximation   is  based   on   the 

use  of  an  ordered   set  of  vector-valued   control  points    for   the  A.'s. 

1 

These  points  represent  a  control  polygon  for  curves  (CP.)  or  a  control 
mesh  (CM.  .)  for  surfaces  that  will  approximate  the  shape  in  the 
geometric  coordinate  system. 


Types  of  B-spline  Representations 

Some  of  the  ways  that  geometrical  applications  of  B-splines  can  be 
categorized  are  by,  1)  dimensionality,  2)  periodicity,  or  3)  spatiality. 
Dimensionally ,  there  are  two  forms  of  B-spline  representations,  the 
integral  and  the  rational. 

The  original  development  of  the  integral  form  was  done  by  Riesenfeld 
(ref  1)  in  1973  and  is  a  natural  extension  of  the  B-spline  approximation 
to  two  or  three  Cartesian  coordinates.   In  his  thesis  he  reviews  the 
mathematical  basis  for  the  B-spline  approximation  and  shows  that  it  is  a 
proper  extension  of  the  Bernstein  polynomial  basis  used  by  Bezier. 
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The  original  development  of  the  rational  form  was  done  by  Versprille 
(ref  2)  in  1975,  in  which  a  homogenous  representation  is  used  to  obtain 
desired  geometric  properties.   The  rational  form  permits  the  exact 
representation  of  conic  sections.  Since  conies  are  the  basis  of  most 
engineering  drawing  and  since  the  integral  form  of  the  B-spline  is  a 
subset  of  the  rational  form  obtained  by  setting  the  homogenous 
coordinate  equal  to  1.0,  the  rational  form  is  better  suited  to  geometric 
applications  in  engineering  design.  The  development  that  follows  will  be 
for  the  rational  form. 

Geometric  periodicity  usually  indicates  whether  the  curve  or  surface  is 
open,  having  distinct  end  points,  or  closed,  and  not  having  distinct 
endpoints.   Open  B-spline  curves  and  surfaces  have  non-periodic  knot 
vectors  with  distinct  beginning  and  ending  sequences  of  multiplicity  K, 
and  the  control  mesh  or  polygon  is  usually  open  with  distinct  starting 
and  ending  points.   Closed  B-spline  curves  and  surfaces  usually  have 
periodic  knot  vectors,  and  the  control  mesh  or  polygon  is  closed  on 
itself. 

Spatiality  refers  to  the  type  of  geometric  entity  that  is  generated. 
B-spline  curves  have  only  one  parametric  variable  and  are  evaluated  by 
summing  the  product  of  the  basis  functions  and  the  active  control 
polygon.   Surfaces  have  two  parametric  variables  and  are  evaluated 
using  a  Cartesian  product  of  basis  functions  in  two  directions  and  the 
elements  of  the  active  control  mesh.   The  integral  form  of  the  B-spline 
curve  is 
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P(t)=sumover(i:j-(K-l),j)  CP .  *  N.   (t) 

X      X  >J\ 

with  i  such  that   t .  <«  t  <  t .  ,  . 

Where  P(t)  and  CP .  are  vectors  in  the  geometric  space. 

The  rational  form  of  the  B-spline  curve  is: 

Q(t)  =  sumover(i:j-(K-l),j)  CPQ.  *  N.  _(t) 

1     1  jK 

with  j  such  that  t .  <=  t  <  t .  .  . 

Where  Q(t)  and  CPQ.  are  vectors  in  the  homogenous  space  and  Q,  (t)  is  the 

i  h 


homogenous  component  of  Q(t).   Then 
P(t)  T  =  |  Q(t)  |  /  Q  (t) 


I E?  I 


B-spline  curves  can  be  applied  in  either  two  or  three  dimensional  space. 


The  integral  form  of  a  B-spline  surface  is 

S(u,v)  =  sumover(i:g-(l-l) ,g)  sumover( j :h-( J-l) ,h)  CM.  .  *  N.  T(u) 

*  N.  _(v) 
J.J 

with  g  such  that  u  <=  u  <  u  ,   and  h  such  that  v,  <=  v  <  v.  , 

g         g+1  h         h+1 

where  S(u,v)  and  CM.  .  are  vectors  in  the  three  dimensional  geometric 

1 » J 

space.   The  rational  form  for  a  surface  is 


R( 


u,v)  ■  sumover(i:g-(l-l) ,g)  sumover( j :h-( J-l) ,h)  CMR.  .  *  N.  T(u) 


*  N.  _(v) 
J.J 


with  g  such  that  u  <=  u  <  u   ,   and   h  such  that  v,  <=  v  <  v,  , ,  . 
6  g         g+1  h         h+1 

Where  R(u,v)  and  CMR.  .  are  vectors  in  the  homogenous  space  and  K(u,v) 

1 .  J  « 

is    the  homogenous   component   of  R(u,v).      Then 
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S(u,v)  t  =  |  R(u,v)  |  /  R^Cu.v) 


T  S(u,v)  t 
I  1.0    J 


Geometric  Properties  of  B-splines 

Both  integral  and  rational  forms  of  the  B-spline  function  exhibit 
similar  geometric  properties.   Since  the  integral  form  is  a  subset  of 
rational  form,  the  development  will  focus  on  the  rational  form. 

Both  curves  and  surfaces  in  the  B-spline  representation  are  invariant 
under  homogenous  transformations,  (rotations,  scaling,  translation,  and 
perspective  transformations)  see  (ref  5). 

There  are  three  different  mechanisms  that  can  be  used  to  control  the 
shape  of  the  B-spline  curves  and  surfaces  corresponding  to  different 
parts  of  the  functional  definition.   They  are  the  order  of  the  curve, 
the  knot  vector  and  the  control  points.   The  relationships  between  these 
mechanisms  is  not  unique,  and  no  intuitive  scheme  is  known. 

The  order  of  the  basis  function  determines  the  minimum  number  of 
control  points  required,  the  overall  shape,  and  the  continuity  (K-2)  of 
the  resulting  curve  or  surface. 

The  knot  vector  provides  a  mechanism  that  can  be  used  to  change  the 
shape  of  curves  and  surfaces,  particularly  by  using  multiplicity  to 
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locally  reduce   the   order  of   the  basis    functions.      The  knot  vector   for  a 
non-periodic   B-spline 

T   :    (t0,tlt    ...    tn+R) 

has  n+K+1  elements  which  satisfy  the  relationship  t.<=t.  ,.  Also, 

1        l+l 

t  =t,=   ...=tL_    -    and   t      .=t   J    =   ...    t    ,_   and    for  uniform  B-splines     with 
0      1  K-l  n+1      n+2  n+K  r 

no  multiplicities    t  ■  n-K+2.   A  common   form  for   the  knot  vector   is 

max 

T  :  (0,0,. ..,0,1,2,  ...  t    ,...,t    ,t   ) 

max      max  max 


K  K 

The  periodic  knot  vector  is 

T:  (t.,t.,...t ,t   .)  with  t.=(i  -  K/2)       . 

0   1     i     n+1        i         mod  n+1 

The  K/2  is  integer  division  and  the  knot  vector  has  no  multiplicity. 


The  points  in  the  control  polygon  or  control  mesh  provide  the  most 
intuitive  control  over  the  shape  of  the  curve  or  surface.   There  is  a 
one  to  one  correspondence  between  the  number  of  control  points  and  the 

number  of  K   order  basis  functions.  Two  interesting  special  cases 
result  from  this  relationship;   one  is  associated  with  multiplicity  of 
control  points  and  the  other  is  associated  with  collinear  control 
points.   If  K  control  points  CP .  are  the  same 

CPi  =  CPi+l=  ••*  =  CPi+(K-l) 
then,  since  the  sum  of  basis  functions  is  always  equal  to  1.0,  the  curve 
must  pass  through  that  control  point.   If  K  control  points  are 
collinear,  then  the  curve  will  be  locally  linear  on  that  span.   If 

CPi'  CPi+l'  '••  '  CPi+(K-l) 
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lie  on  some  line,  then  the  curve  spanned  by  those  control  points  will  be 
linear.   These  special  cases  are  the  result  of  two  properties  of  the 
B-spline  function.   The  first,  is  variation  diminishing,  and  second,  it 
has  local  support.   These  two  properties  combined  give  rise  to  the 
property  described  later  as  the  convex  hull. 

The  geometric  relationships  between  these  control  mechanisms  are  not 

unique  In  Rogers  and  Adams  (ref  5)  a  method  is  described  that  attempts 

to  directly  relate  repeated  control  points  to  multiplicities  in  the  knot 

vector;  in  Tiller's  article  (ref  6),  there  are  multiplicities  in  the 

knot  vector  without  multiple  control  points.   In  general  the  knot  vector 

must  be  specified  such  that  the  N.  _(t)  are  correctly  defined  with 

ljK. 

respect  to  the  definition  of  the  B-spline  approximation. 


B-spline  Curves 

Geometric  properties  of  interest  in  the  application  of  curves  to 
engineering  design  are:  1)  local  support,  2)  convex  hull,  3)  local 
linearity,  and  4)  numerical  evaluations.   Local  support  means  that  there 
will  be  at  most  K  control  points  required  to  evaluate  any  point  t  on  the 
curve.   These  control  points  are  referred  to  as  the  active  control 
points  or  the  supporting  polygon.   The  convex  hull  property  states  that 
any  point  on  the  curve  will  lie  within  the  convex  hull  of  the  supporting 
polyhedron.   Figure  3-1  illustrates  the  convex  hull  property  for  a 

planar  curve.   Local  linearity  implies  that  within  any  K   order  curve 
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there  can  be  locally  linear  portions  if  there  are  K  collinear  control 
points.  Figure  3-2  illustrates  local  linearity. 

Two  vectors  associated  with  geometric  properties  of  interest  are:  the 
position  vector  and  the  tangent  vector.   The  position  vector  is  obtained 
for  the  rational  B-spline  by  dividing  by  the  homogenous  component. 

The  formula  for  the  evaluation  of  points  on  a  B-spline  curve  was  given 

at  the  beginning  of  this  section.   The  evaluation  of  the  tangent  vector 

to  the  B-spline  curve  is  the  vector  extension  of  the  derivative  of  the 

B-spline  approximation  derivative,  with  the  derivative  constants 

becoming  a  set  of  derivative  control  points  evaluated  in  vector  form  for 

d=l. 

To  evaluate  the  rational  B-spline  tangent  vector: 

first  evaluate  the  parametric  derivatives  in  homogenous  form 

Q(1)(t)  =  sumover(i:j-(K-l)+l,j)  CPQ(^  *  N£  K-1(t) 

j  s.t.  t.  <=  t  <t.  . . 

Then  the  geometric  tangent  vectors  are 

"  -  |  Qh(t)*Q(1)(t)  -  Q(t)*Q^1}(t)  |  /  Qh2(t). 


i  ■:::<■'  \ 


The  B-spline  approximation  is  defined  as  being  left  continuous  for  the 
interval  t.<=  t  <  t.  ..   For  open  curves,  as  the  parameter  approaches 

the  maximum  value  of  the  knot  vector  t    ,  all  the  basis  functions 

max 

approach  zero  except  N  ^(t)  which  approaches  1.0.   The  curve  is  equal 

n  ,K. 

to  the  end  point. 
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B-spline  Surfaces 

The  geometric  properties  of  interest  in  application  of  B-spline  surfaces 
to  engineering  design  are  similar  to  those  for  curves.   These  are:  local 
support,  convex  hull,  locally  planar  surfaces  and  numerical  evaluations. 
For  a  surface  there  are  two  parametric  variables  u  and  v.   Each  of  the 
parameters  is  associated  with  a  basis  function  which  has  an  order,  a 
knot  vector,  and  a  control  mesh  instead  of  a  polygon. 

Local  support  for  a  surface  means  that  for  any  point  (u,v)  on  the 
surface  there  will  be  a  control  mesh  that  will  be  at  most  I  by  J  points 
needed  to  evaluate  it  and  its  derivatives.   I  is  the  order  in  the  u 
direction  and  J  is  the  order  in  the  v  direction.   The  mesh  will  be 
called  the  active  mesh  or  supporting  polyhedron. 

The  convex  hull  for  a  surface  is  defined  as  the  space  bounded  by  the 
supporting  polyhedron,  and  a  point  on  the  surface  will  lie  within  that 
space.  Figure  3-4  shows  the  convex  hull  property  for  a  surface. 
Portions  of  a  surface  can  be  made  to  be  locally  linear  by  making  I 
control  points,  in  the  u  direction  or  J  control  points,  in  the  v 
direction  collinear.  Portions  of  the  surface  can  be  made  to  be  locally 
planar  by  making  an  I*J  array  of  control  points  coplanar. 

The  two  properties  required  for  surface  applications  are  the  position 
vector  and  the  unit  normal  to  the  surface  at  the  point.   Points  on  the 
B-spline  surface  are  evaluated  using  a  Cartesian  product  of  the  active 
control  mesh  and  the  appropriate  basis  functions  in  each  direction. 
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This  was   shown  at  the  beginning  of  this  section  for  both  the  rational 
and  integral  form. 

Unit  normals  to  the  surface  are  evaluated  by  taking  the  cross  product  of 
the  geometric  tangent  vectors  in  the  u  and  v  directions.  The  tangent 
vectors  can  be  evaluated  using  the  parametric  partial  derivatives  of 

R(u,v).   R     is  the  partial  derivative  of  R(u,v)  with  respect  to  u, 

and  R   '    is  the  partial  derivative  of  R(u,v)  with  respect  to  v. 

R(u,v),  R   '   ,R      are  homogenous  vectors.   To  evaluate  the  partials 

the  homogeneous  derivative  control  meshes,  CMR.  *.    and  CMR.  '.   ,  are 

used.   The  tangent  vectors,  S   '   and   S        and  the  normal  vector, 
NV,  are  geometric  vectors.   The  procedure  to  evaluate  the  unit  normal  to 
the  surface  is: 
1.    Evaluate  the  parametric  partial  derivatives  of  R(u,v). 

a.  evaluate  the  derivative  control  points  in  each  direction 

CMR.!v)  =  CMR. 

for  the  partial  of  R(u,v)  with  respect  to  u 

CMR^!0  =  CMRU!V)-  CMR^V-  /  (u.  T  ,  -  u.) 
i,j         i,j        l-l  ,j      l+I-l     l 

(i:g-(I-l)+l,g)  (j:h-(J-l),h) 

and    for   the   partial   of   R(u,v)   with   respect    to  v 

CMR<V!°    =  CMRU!V)-  CMRU!VJ    /    (v.    f   .    -  v.) 
x,j  i,j  i>J~l  i+J-1  1 

(i:g-(I-l),g)    (j:h-(J-l)+l,h). 
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b.  evaluate  the  Cartesian  product  of  derivative  control  points 
and  the  reduced  order  basis  functions  for  each  partial. 

R   "    ■  8umover( i:g-(l-l)+l ,g)sumover( j :h-( J-l) ,h) 

CMR^U!1)  *  N.  _  ,(u)  *  N.  T(v) 
i,j      1,1-1       j, J 


R      ■  sumover(i:g-(l-l) ,g)sumover( j :h-( J-l)+l ,h) 

CMR^U!U  *  N.  _  .(u)  *  N.  _(v) 
1  |  j       1  >  l— l        2  » •* 


2.   Calculate  the  geometric  tangent  vectors  S     and  S 


j £■»  j.  t«<«».H -..«£«»!/.,,' 


(v:l)  t 


S 

0.0 


(v:l) 


.j,^».v,.^j,h 


3.   Evaluate  the  cross  product  to  obtain  the  geometric  normal  to  the 
surface. 


•■  l 


k 


S(u:1)  X  S(v:1)  ■■■   \   c(u:l)   _(u:l)   c(u:l) 

D  D  S 

s(v:l)   s(v:l)   s(v:l) 
+  x       y       z 

■  (NV  ,  NV  ,  NV  ) 
x    y    z 

4.    Normalize  the  normal  vector  to  get  the  unit  normal  to  the  surface 

at  R(u,v). 

Unit  Normal  Vector  =  NV  /  SQRT(NV2  +  NV2  +  NV2). 

x     y     z 
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For  open  surfaces  two  special  boundary  cases  exist.   First,  when  either 
of  the  parameters  is  equal  to  the  maximum  knot  value  in  that  direction, 
the  supporting  mesh  degenerates  to  a  supporting  polygon  along  the  edge 
in  the  other  direction.   Figure  3-4  shows  this  case.   Second,  when  both 
parameters  are  equal  to  the  maximum  knot  values,  the  supporting  mesh 
degenerates  to  the  corner  control  point. 

This  completes  the  discussion  of  the  geometric   application  of  the 
B-spline  approximation.   The  next  section  will  describe  the  program  that 
was  used  to  display  and  modify  B-spline  surfaces. 
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IV.  A  PROGRAM  TO  GENERATE,  MODIFY,  AND  DISPLAY  3-D  B-SPLINE  SURFACES 


The  Implemented  Structure 

The  program  BSPLINE  is  an  interactive,  menu  guided  program  that  can 
generate,  display,  and  modify  integral  and  rational  non-periodic 
B-spline  surfaces.   There  are  three  sets  of  arrays  used  by  the  program, 
database  arrays,  the  display  arrays,  and  the  display  file  arrays.   The 
database  arrays  contain  the  data   from  the  data  file  which  includes:  1) 
the  order  of  the  basis  functions,  2)  the  number  of  the  control  points  in 
each  direction,  3)  the  knot  vectors,  4)  the  dimension  of  control 
points  and  5)  the  coordinate  vectors  of  the  control  points.  The  display 
arrays  contain  information  used  to  display  the  streamlines  and  normals 
to  the  surface  in  the  GENERATE  BSPLINE  module  of  the  program  described 
later.  The  display  file  arrays  contain:  1)  the  index  counters  needed  to 
display  a  display  file,  2)  the  geometric  location  of  the  control  points, 
and  3)  the  information  needed  to  regenerate  the  streamlines  and  normals 
to  the  surface.  These  are  used  by  the  module  DISPLAY  FILE  MANAGER  (to  be 
described  later)  to  save  and  read  data  from  external  storage. 

The  overall  structure  of  the  program  is  shown  in  figure  4-1.   The  main 
program  controls  access  to  the  different  modules  in  the  program.   The 
module  GENERATE  BSPLINE  is  used  to  generate  the  streamlines  and  normals 
and  to  orient  the  control  mesh.   The  submodule  MODIFY  BSPLINE  is  used  to 
locally  modify  B-spline  control  mesh.   The  DISPLAY  FILE  MANAGER  is  used 
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to  save  and  display  files.   Appendix  B  contains  an  index  of  the 
subroutines  used  by  the  program  and  a  subroutine  call  tree  organized 
around  the  program  structure. 


Ha in  Program  Module 

BSPLINE  MAIN  serves  as  a  control  structure  to  enter  the  modules  that 
have  been  implemented  and  to  enter  and  exit  the  main  program.   The  menu 
and  presentation  window  are  shown  in  Figure  4-2.   The  commands  that  are 
available  from  the  main  program  are 

GENERATE  BSPLINE 

DISPLAY  FILE  MANAGER 

PRINT  FILE 

HELP 

GENERATE  BSPLINE  is  a  module  that  is  used  to  orient  the  control  mesh  for 
viewing  and  to  generate  the  the  display  arrays.  A  more  complete 
description  will  be  given  in  a  following  section. 

DISPLAY  FILE  MANAGER  is  used  to  create,  save,  and  display  the  display 
files.  A  more  complete  description  of  this  will  also  be  given  later. 

PRINT  FILE  is  an  interactive  segment  that  permits  the  user  to  look  at 
the  current  program,  database,  and  display  arrays.   This  command  is 
available  from  anywhere  in  the  program  and  is  on  all  menus.   The  primary 
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function  is  to  provide  dynamic  feedback  from  the  program  during  the 
development  of  the  different  modules.   The  information  displayed  can 
also  be  written  to  system  files.   A  different  file  can  be  created  each 
time  the  command  is  executed,  \ 

HELP  provides  information  for  any  command  in  the  current  menu  by  typing 
<HP>  <menu  item>.  If  no  menu  item  is  listed,  the  structure  of  the  help 
command  is  shown.   This  command  was  not  implemented  in  the  present 
version  of  the  program. 

These  are  the  commands  available  from  the  main  menu  of  the  program. 


Module  to  Generate  B-spline  Surfaces 

GENERATE  BSPLINE  is  an  independent  module  that  inputs  data  files, 
orients  the  control  mesh  for  viewing,  and  generates  the  arrays  used  to 
display  the  surface.   The  menu  and  presentation  window  are  shown  in 
figure  4-3.  The  commands  available  from  this  menu  are 

READ  FILE 

DRAW  MESH 

DRAW  NORMALS 

DRAW  STREAMLINES 

ISOMETRIC 

ROTATE 

SCALE  DATA 
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STREAMLINES  AND  NORMALS 

MODIFY  BSPLINE 

PRINT  FILE 

ERASE 

RESET 

SAVE  FILE 

QUIT 


READ  FILE  and  SAVE  FILE  are  used  to  read  and  save  the  database  data 
from  external  storage.   The  user  specifies  the  <filename>  and  <filetype> 
of  the  external  file.   Examples  of  data  files  are  given  in  Appendix  C 
and  include  the  orders,  the  knot  vectors,  the  number  of  vertices  in  each 
direction,  the  dimension,  and  the  coordinates  of  the  control  mesh. 

STREAMLINES  AND  NORMALS  is  used  to  generate  the  streamlines  and  normals 
in  the  display  arrays  used  to  display  the  surface.   The  user  is  queried 
for  the  number  of  streamlines  and  points  to  generate  per  interval  on  the 
surface.   The  normals  are  generated  at  intersections  of  streamlines. 
Current  values  for  the  streamlines  and  normals  must  exist  in  order  to 
use  the  drawing  commands  and  to  enter  the  MODIFY  BSPLINE  module  of  the 
program. 

DRAW  MESH,  DRAW  NORMALS,  and  DRAW  STREAMLINES  are  used  to  draw  the 
respective  parts  of  the  database  and  the  display  arrays  if  they  exist. 
DRAW  MESH  draws  the  control  mesh  shown  in  figure  4-4.   DRAW  NORMALS 
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draws  the  unit  normals  at  the  intersections  of  the  streamlines. 
Figure  4-5  shows  the  normals.   DRAW  STREAMLINES  draws  the  streamlines 
specified  in  STREAMLINES  AND  NORMALS.   Figure  4-6  shows  an  example  of 
the  streamlines.   Any  attempt  to  draw  something  that  does  not  exist 
generates  a  warning  from  the  program. 

ISOMETRIC,  ROTATE,  and  SCALE  DATA  are  used  to  modify  the  control  mesh  in 
the  database  for  viewing.  ISOMETRIC  generates  an  isometric  view  by 
rotating  the  control  mesh  by  45  degrees  about  the  y  axis  followed  by 
35.246  degrees  about  the  x  axis.  An  example  of  an  isometric  view  is 
shown  in  figure  4-7.  ROTATE  rotates  the  control  mesh  in  the  database 
about  an  <axis>  by  an  angle  specified  by  the  user.  Multiple  angles  can 
be  specified.  The  rotations  are  performed  in  the  sequence  given. 
Figure  4-8  shows  a  data  set  that  has  been  rotated  about  x  by  60  degrees, 
about  y  by  45  degrees  and  about  z  by  10  degrees.   SCALE  DATA  is  used  to 
scale  the  control  mesh  in  the  database  by  a  factor  <scl>  specified  by 
the  user.   The  factor  must  be  greater  than  .001.  Figure  4-9  shows  that 
both  the  old  data  and  the  scaled  data  are  displayed.   Each  time  these 
commands  are  used,  the  user  is  queried  whether  or  not  to  replace  the 
current  data  with  the  modified  data.  If  the  modified  data  file  is  used 
the  streamlines  and  normals  must  be  regenerated. 

MODIFY  BSPLINE  is  a  module  used  to  locally  modify  the  control  mesh  and 
can  only  be  entered  if  current  streamlines  and  normals  exist.   This 
module  will  be  described  in  the  next  section. 
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PRINT  FILE  was  described  in  the  discussion  of  the  main  program. 

ERASE,  RESET,  and  QUIT  are  local  utilities  used  to  erase  the  screen, 
reset  the  database,  and  quit  this  module.   Only  one  data  file  can  exist 
in  the  database;  therefore,  the  database  must  be  reset  before  a  new  data 
file  can  be  read. 


Display  File  Module 

DISPLAY  FILE  MANAGER  is  an  independent  module  that  is  used  for  the 
display  of  the  control  mesh,  streamlines  and  normals.   It  both  reads 
and  saves  display  files  from  memory  or  creates  them  from  the  display 
arrays  used  in  GENERATE  BSPLINE  if  they  exist.   The  menu  and 
presentation  window  are  shown  in  figure  4-10.   The  commands  available 
from  this  menu  are: 

CREATE  FILE 

READ  FILE 

DRAW  MESH 

DRAW  STREAMLINES 

DRAW  NORMALS 

SCREEN  TYPE 

SCALE  OUTPUT 

PRINT  FILE 

ERASE 

RESET 
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SAVE  FILE 
QUIT 

CREATE  FILE  creates  a  display  file  by  copying  the  control  mesh,  the 
display  arrays,  and  the  counters  needed  to  display  them  into  the  display 
file  arrays  from  GENERATE  BSPLINE  arrays.   These  arrays  can  then  be 
saved  on  external  media. 

READ  FILE  and  SAVE  FILE  are  used  to  read  and  save  display  files  to  and 
from  external  storage.   The  data  file  can  not  be  constructed  from  the 
data  in  the  display  file. 

DRAW  MESH,  DRAW  NORMALS,  DRAW  STREAMLINES  are  similar  to  the  commands 
given  in  GENERATE  BSPLINE  except  they  use  the  display  file  arrays. 

SCREEN  TYPE  permits  the  user  to  specify  the  presentation  window  as 
either  the  standard  window  with  menu  or  a  full  screen  without  menu. 
Figure  4-11  shows  a  full  screen  display. 

SCALE  OUTPUT  scales  the  size  of  the  output  window  when  using  the  full 
screen  mode.  Figure  4-12  shows  a  20%  reduction  in  the  size  of  the 
window. 

PRINT  FILE  is  the  same  as  in  the  main  program. 
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ERASE,  RESET,  AND  QUIT  are  utilities  used  to  clear  the  screen, 
initialize  the  display  arrays,  and  quit  this  module. 

These  are  the  commands  available  in  the  display  module. 

General  Information 

The  program  consists  of  approximately  7,000  lines  of  commented  FORTRAN 
77  source  code.   It  was  written  on  the  Kansas  State  University  Computing 
Center's  MAS  6630  mainframe  using  a  subset  of  PLOT  10  IGL  graphics 
subroutines.   The  output  devices  were  a  Selenar  Hirez  100  graphics 
terminal  setup  to  dump  to  a  Hewlett/Packard  7470A  plotter. 

In  this  chapter  the  operation  of  the  MAIN  program,  GENERATE  BSPLINE,  and 
DISPLAY  FILE  MANAGER  modules  have  been  explained  and  illustrated.  The 
next  chapter  explains  the  operation  of  the  MODIFY  BSPLINE  module. 
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BSPLINE  MAIN 
MODULE 

GENERATE  BSPLINE  <GNB >+++++++ 

DISPLAY  FILE  MANAGER  <DFM>++ 

PRINT  FILE  <PF> 

HELP  <HP> 

QUIT  PROGRAM  <QT> 


I 
I 
I 
I 
V 


GENERATE  BSPLINE 
MODULE 

READ  FILE  <RF> 

DRAW  MESH  <DM> 

DRAW  NORMALS  <DN> 

DRAW  STREAMLINES  <DS> 

ISOMETRIC  <IS> 

ROTATE  <RT> 

SCALE  DATA  <SD> 

STREAMLINES  AND 

NORMALS  <SN> 

MODIFY  BSPLINE  <MDB>++++ 


DISPLAY  FILE 
MANAGER  MODULE 

CREATE  FILE  <CR> 
READ  FILE  <RF> 
DRAW  MESH  <DM> 
DRAW 

STREAMLINE  <DS> 
SCREEN  TYPE  <ST> 
PRINT  FILE  <PF> 
ERASE  <ER> 
RESET  <RS> 
SAVE  FILE  <SV> 


PRINT  FILE   <PF> 

|           QUIT   <QT> 

ERASE   <ER> 

V 

RESET   <RS> 

MODIFY  BSPLINE 

SAVE  FILE   <SF> 

MODULE 

QUIT   <QT> 

MOVING  POINTER  <MP> 

DRAW  ACTIVE  MESH  <DA> 

DRAW  MESH  <DM> 

DRAW   STREAMLINES   <DS> 

MOVE  PIJ"S   <MV> 

UPDATE   <UP> 

PRINT  FILE   <PF> 

ERASE   <ER> 

QUIT   <QT> 

Figure  4-1  Program  structure 
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V.  INTERACTIVE  MODIFICATION  OF  B-SPLINE  SURFACES 


MODIFY  BSPLINE  is  a  module  that  can  be  accessed  from  GENERATE  BSPLINE 
when  current  streamlines  and  normals  exist.  The  module  is  used  to 
locally  modify  the  surface  by  moving  the  points  of  the  control  mesh  in 
space.  The  menu  and  presentation  window  are  shown  in  figure  5-1.  The 
commands  available  from  this  menu  are: 


MOVING  POINTER 

DRAW  ACTIVE  MESH 

DRAW  MESH 

DRAW  STREAMLINES 

MOVE  PIJ"S 

UPDATE 

PRINT  FILE 

ERASE 

QUIT 


Relevant  Geometric  Properties 

The  geometric  properties  that  are  used  in  this  module  are  the  local 
support,  the  convex  hull,  and  the  existence  of  a  u,v  surface  coordinate 
system. 
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Local  support  guarantees  that  there  will  be  a  finite  mesh  of  control 
points  associated  with  the  generation  of  any  point  on  the  surface.  The 
active  control  mesh  will  be  an  array  of  control  points  no  larger  than 
the  order  of  the  surface  in  the  u  direction  by  the  order  of  the  surface 
in  the  v  direction. 

The  convex  hull  property  of  the  B-spline  surface  guarantees  that  any 
point  on  the  surface  will  be  in  the  space  of  the  convex  hull  of  the 
polyhedron,  defined  by  the  active  control  mesh. 

A  u,v  coordinate  system  results  from  the  use  of  the  parametric 
representation  of  the  surface.  The  parameters  map  the  three-dimensional 
x,  y,  z  geometric  space  onto  the  two-dimensions  of  the  parameters  u  and 
v.  This  mapping  permits  the  effective  use  of  two  dimensional  input 
devices  such  as  thumb  wheels,  cursor  keys,  or  a  mouse  to  move  a  place 
marker  on  the  surface. 

The  combination  of  these  properties  led  to  the  development  of  an 
effective  means  of  positioning  a  cursor  on  the  surface.   Associated  with 
that  position  is  an  active  control  mesh,  with  which  to  modify  a 
that  portion  of  the  curve. 
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Critical  Conmands 

MOVING  POINTER  is  used  to  move  a  surface  pointer  on  the  u,v  surface  and 
to  dynamically  display  the  active  control  mesh  associated  with  the 
point.  Pointer  position  is  controlled  using  four  cursor  keys,  +u,  -u, 
+v ,  -v . 

The  one-to-one  relationship  between  a  point  on  the  surface  and  the 
active  control  mesh  required  to  generate  it  permits  the  dynamic  display 
of  the  active  control  mesh.   Shown  in  Figure  5-2.  Both  the  parametric 
and  the  Cartesian  location  are  displayed  and  updated  with  each  movement 
of  the  pointer.  The  active  control  mesh  is  displayed  when  DRAW  ACTIVE 
MESH  is  on.  The  purpose  of  MOVING  POINTER  is  to  position  the  pointer  on 
the  surface  and  establish  the  mesh  that  can  be  modified  in  MOVE  PIJ"S. 

MOVE  PIJ"S  is  used  to  change  the  location  in  space  of  the  elements  of 
the  active  control  mesh.  The  active  control  mesh  is  determined  by  the 
location  of  surface  pointer  in  MOVING  POINTER.   An  interactive  sequence 
precedes  the  modification  of  control  points  where 

1.   The  user  chooses  whether  to  display  the  entire  control  mesh  or 
zoom  the  active  control  mesh. 

If  the  user  chooses  to  zoom  the  active  control  mesh  then 
only  the  active  control  mesh  and  the  adjacent  control 
points  are  displayed. 
(Figures  5-4  and  5-5  illustrate  this  difference.) 
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2.  The  user  sets  the  mover  sensitivity,  a  control  which  adjusts 
the  incremental  amount  the  points  are  moved  for  each  stroke  of 
the  control  keys.  The  default  is  .01  of  the  maximum  coordinate 
difference  in  the  x,  y,  z,  directions  that  the  mesh  occupies. 

3.  The  labels  are  listed,  and  the  user  chooses  one  of  the  elements 
of  the  active  control  mesh  to  modify  by  specifying  its  I, J  mesh 
index. 

Six  control  keys  are  used  to  change  the  position  of  the  selected  point: 
+x,  -x,  +y,  -y,  +z ,  -z.  Multiple  points  can  be  changed,  or  the  user  may 
exit  the  command. 

UPDATE  is  used  to  recalculate  the  streamlines  and  normal  after   the 
control  mesh  is  modified.  The  point  and  interval  counters  are  the  same 
as  those  used  in  STREAMLINES  AND  NORMALS  in  GENERATE  BSPLINE  module. 
The  display  arrays   must  be  updated  after  any  modification  of  the 
control  mesh  to  be  able  to  display  the  streamlines. 


Other  Commands 

DRAW  ACTIVE  MESH  is  used  to  activate  the  dynamic  drawing  of  the  active 
mesh  in  MOVING  POINTER. 

DRAW  MESH  and   DRAW  STREAMLINES  are  used  to  draw  the  current  mesh  and 

s  t reaml inet . 
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PRINT  FILE   is  the  same  as  in  the  main  program. 

ERASE  is  a  utility  to  erase  the  screen. 

QUIT  quits  this  submodule  and  returns  to  GENERATE  BSPLINE. 

This  completes  the  description  of  the  module  to  locally  modify  the 
B-spline  control  mesh.  Next  is  a  review  of  what  has  been  done  with 
possible  recommendation  for  future  work. 
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VI.  CONCLUSIONS 


Accomplishments  of  this  present  work 

There  were  three  main  accomplishments  in  this  work:  1)  the  researching 
of  the  B-spline  surface  representations,  2)  the  development  of  an 
interactive  program  to  generate  and  display  B-spline  surfaces,  and  3) 
the  implementation  of  a  method  for  local  modification  of  B-spline 
surfaces. 

The  research  involved  an  extensive  review  of  both  current  and  historic 
literature  about  the  application  of  the  B-spline  approximation  to  the 
computer  generation  of  surface  models.   Using  the  information  aquired  in 
the  literature  search,  an  interactive  program  was  written  to  generate 
and  display  both  integral  and  rational  B-spline  models.  Finally,  to 
extend  the  usefulness  of  the  B-spline  representation,  a  method  of 
locally  modifying  B-spline  surfaces  was  developed.   The  basis  for  the 
development  was  the  local  support  and  convex  hull  properties.   The 
method  utilized  the  inherent  two-dimensional  characteristics  of  the 
parametric  representation  to  move  a  location  pointer  on  the  surface. 
Using  the  location  of  the  parametric  pointer  to  specify  a  finite  set  of 
control  points,  an  intuitive  method  of  modifying  the  shape  of  the  curve 
was  developed. 
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Possible  Extensions  of  this  Work 

Some  directions  that  could  be  taken  in  extending  this  work  are:  further 
development  of  the  program,  use  of  the  present  program  as  an  educational 
tool,  and  the  development  of  applications  based  on  the  surface 
representation. 

Other  modules  that  would  be  useful  are 

1.  an  interactive  database  builder  that  would  permit  the 
construction  of  original  B-spline  surfaces 

2.  a  module  to  permit  the  modeling  of  multiple  surfaces  at  the 
same  time 

3.  a  module  for  large  scale  reshaping  of  the  control  mesh 

4.  the  development  of  a  more  sophisticated  display  module  to 
include  hidden  surface  removal  and  shadowing  algorithms  to  be 
used  when  displaying  surfaces. 

The  present  program  and  thesis  could  be  used  as  an  educational  tool  to 
introduce  students  to  the  characteristic  of  the  B-spline  surface  and  to 
the  algorithms  required  to  generate  them.   The  generate  module  could  be 
used  in  another  program  to  model  approximate  multi-dimensional  data 
sets. 

Various  modules  developed  in  this  work  could  be  used  individually  or 
together  in  graphical  applications  for  mechanical  part  design,  analysis, 
or  manufacture. 
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The  B-spline  representation  for  graphical  entities  is  a  powerful 
mathematical  tool  in  the  development  of  engineering  analysis  and  design. 
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APPENDIX  A 
B— spline  Basis  Function  Evaluations 
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B-SPLIHE  FCTPEMEKEAL  EVALC&TICIg 
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therefore  N.  ,(t):  not  equal  to  0.0 
i,l  ^ 


H.  .(t)=  1. 
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If  K^  then  sxmover(i:j-(K-l),j)=sumover(i:j-l,j) 

j  s.t.  t.  <=  t  <  t.  , 
J  J+l 
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Ifote:   i=i-l  in  the  general  definition 
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Note:  i=i  in  tte  general  definition 
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n(i-3)  Ki-2)^ 


H.  .  ,(t)  =  t.       -  t 
1-3,4  l+l 


'i+l  "  Ci-2 


'irt-* 


'i+1  "  Ci-1 


N.  ,(t) 


+  + 
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Nate:  i=i-2  in  the  general  definition 


Ni-2,4(t)  =  fc  -  Hi-l) 


'(1-2) +3  "  '(1-2) 


N(i-2),3  *  t(i-2)j4  "  fc 


N(i-2)+l,3(t) 


'(1-2)44 -'(1-2)+! 


I(i-3):H(i-2),3 


II(i-2):H(i-l)^ 


"t-t(H)        N(i-2),3  +  ti12-i_     Ni-l,3(t) 


'i+1  '  '(1-2) 


Ci+2  "  Vl 


l(i-2)  H(i-2),3 

x-2,4 


C  -  '(!-» 

.  'i+1  "  'i-1 

t.  -   -  t    N.  .(t) 
l+l               1,1 

Ci+1  "  C(i-2! 

• 

n(i-2) 

tfi-4),3a 
'"'(i-l) 

h+i " '  Ni,i(t) " 

e«  -  ci-l 

• 

.  fci+l  ^(i-l)  . 

4- 

n(i-2) 

'i+2  -  t 

'«  -  Vi 


Ki-1)3> 

fci+2-  t 
^2"  fci 


t  -  t.       N.  .(t) 


+  + 


WafaMfeMaMHHfe^^ 


-77- 


Note:  i=i-l  in  the  general  definition 


N.  .   .(t)  =  t  -  t,.  n 
l-l  ,4  (i-l) 


N(i-D,3(t)+t(i^4-l 


N(i-D+l,3(t) 


'(i-D+a  "  c(i-i) 


H i-l)-t4  "  '(i-D+l 


Ki-4):H(i-l),3 


H(i-1):H  i,3 


-  t  -  t,.  .»        ft,       .  _(t)  +  t.      -  t    N.     (t) 
(i-l)  (i-l),3  i+3  i,3 


Ci+2  "  C(i-1) 


t.,  -  t. 
i+3        l 


Ki-l)  tfi-D,3a 


Wt}  =  fc  -  Vd 


Ci+2  "  '(i-l) 


t  -  t 


(i-l) 


Ci+1  "'(i-l) 


fci+l  -  '    Ni,l(t) 
fcM  -  Ci 


+  + 


Ki-4) 

t  -  t 


(i-l) 


fci+2  '  '(i-l) 


*i-l),3b 


t  -  t.      N.  .(t) 

t.  .-  t. 
l+l       l 


n(i-i)       nu 


'i+3  -  t 

fci+3  ;  h 


t  -  t. 

1 

t..,     t. 

i+2  -    i 


t  -  t.       N.  .(t) 


+  + 


*MMaHHHfefcfafcfc^^ 
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Note:  i=i  in  the  general  definition 


I(i):Hi,3 


H(i):tfi+1),3 


N.  ,(t)=t-t.         N.  ,(t)  +  t.   .   -  t        N.  .   _(t) 
i,4  i         i,3  i+4 i+l,3 


'i+3  "  h 


Ci+4  "  ti+l 


n(i) 

Ci+4  "  fci+l 


N(i+1),3 

'  '  fc(i+l) 


N(i+l),2(t)  + 


V+D+S  "  '(i+l) 


fc(i+l)+3  "  t 


N(i+l)+l(t) 


t(i+l)+3  "  ^i+D+l 


I 
V 


II(i) 

Ci+4  "  'i+1 


Ki+D^ 
I(i+l)a:H(i+l),2 


Ki+D^ 
n(i+l)b:H(i+2),2 


t  '  fc(i+l)        N(i+l),2(t)+tij4-i_    Ni+2,2(t) 


Ci+4  "  '(i+l) 


fci+4  "  fci+2 


n(i) 

Ci+4  "  fci+l 


I(i+l)a 

t  -  t 


(i+D 


'i+4  "  c(i+l) 


«i+l),2 
fc  "  C(i+D 


C(i+1)+1  "  '(i+D 


N(i+l),l(t)  + 


t(i+l)+2  "  fc 


N(i+l)+l,l(t) 


t(i+l)+2  "  '(i+D+l 


i  ii  ii  1 1 1 1 1 1 1 1 1 1 1  ti  h  1 1 1 1 1 1 1 1 1 1 1 1 1  ii  nil  1 1 1  ii  1 1 1 1 1  in  1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

H(i+1),3 
H(i)  I(i+l)a 

t  -  t, 


fci+4  "  fci+l 


'(i+D 


Ci+4  "  ^i+l) 


'-'(i+l)       N(i+l),l(t)  + 
Ci+2  "  fc(i+l) 


h+3  -  fc        Ni+2,l(t) 
fci+3  "  Ci+2 
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n(i) 

Ci+4  "  ti+l 


n(i+l)b:Ki+2)^ 

VJ        Ni+2,2(t) 
fci+4  "  ti+2 


II(i) 


*i+l),3 

n(i+i)b 


H(i+2),2 


t.  .  -  t 
i+4 

t.  .  -  t 

i+4 

t'ttm        Wi* 

ti+4  "  fci+l      fci+4  "  ti+2 

,  t(i+2)+l  "  C(i+2) 

C(i+2)+2  "  C             N(i+2)+l,l(t) 

t(i+2)+2  "  C(i+2)+l 

1 
V 
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fci+4  "  ti+l 


n(i+i)b 

ti+4  "  ti+2 


N(i+2),2 

'"'(I*)        N(i+2),1* 
Ci+3  "  t(i+2) 


h*Zl       Ni+3,l(t) 
ti+4  ~  Ci+3 
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I(i)  N  i,3 


N.   ,(t)  -  t  -  t. 
1.4  1 


t..,  -  t. 

i+3        l 


t  -  t. 


t.  -      t. 
1+2  -    l 


t  -  t.       N.     (t) 
i         1,1 

t.  .-  t. 
1+-1       l 


JHcfcfcfcMHHHHMri^^ 
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SQftfAK  (F  BASIS  FWCTICRS 

K=l 

H.   _(t)=  1. 
1*1 

K=2 

"H^(t)  "  h-A  -  '    Hi,l(t) 


B.  ,(t)  =  t  -  t.      N.  ,(t) 
1,2  l        1.1 


K=3 


n(i-2) 

N.       _(t)  =  t.  .   -  t 
1-2,3  l+l 


Ci+1  "  Vl 


N(i-1),2 

t.  .   -  t    N.     (t) 
l+l 1,1 

t.  .   -  t. 
l+l         l 


Ki-4) 
Vl.3(t)  =  t  -  *(£-!) 


'i+l  "*(i-D 


fci+l  -  fc   Ni,l(t) 
fci+l  "  h 


n(i-i)   A  h  i.,2 

t  -  t.       N.  .(t) 
i         1.1 


+  Ci+2-  C 


1+2       l 


'i+r  ci 


+ 


H.  _(t)       t  -  t. 
1.3        _  i 


I(i)  HU 


t.  .       t. 
1+2  -    l 


t  -  t.       N.  .(t) 
i         1,1 

'i+r  ci 


MHHHHfcfcMdfefcfc^^ 
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SDMffif  OF  BASIS  FUNdTOB  (continued ) 
K=4 


H(i-3) 

1-3,4  l+l 


t.  .  -  t.  0 
l+l         1-2 


Ki-2)3 

Vi ; fc 

fcw  -  Ci-1 


tj+l  - fc  Ni,l(t) 


i-2,4 


Ki-2) 
(t)  =  t  -  t 


(i-2) 


fci+l  "  C(i-2) 


tfi-2),3 

t.  .  -  t 
l+l 

t.  .  -  t.  . 
l+l        l-l 


'i+i  -  c  Ni,i(t) 


+  + 


+  + 


nU-2) 

Ci+2  -  C 
h+2  "  Ci-1 


H(i-l),3a 

C  -  C(i-1) 


'w  - C  Ni,i(t) 


n(i-2) 
'i+2  -  h-1 


«i-4)^b 
»«"  'i 


t  -  t.       N.     (t) 
l        1,1 

ci+r  ci 


+  + 


Ki-1) 

H.  ,   .(t)  =  t  -  t(.  n 
i-l,4  (l-l) 


fci+2  "  C(i-1) 


*i-4),3a 

ci+i  _t(i-i) 


fci+i  -  fc  Ni,i(t) 
fci+i  -  \ 


+  + 


i(i-4) 

t  ~  '(i-l) 

h+2 "  '(w) 


Ki-l),3b 


t  -  t.   N.  .(t) 

i    1,1 


n(i-4) 

're  - ■= 

'i+3  "  ei 


Hi,3 
t  -  t. 

^+7^1 


t  -  t.   N.  .(t) 
i    1,1 

ei+r  'i 


+  + 


I(i) 

N.  .(t)  =  t  -  t. 
i,4     i 


'i*  "  Ci 


Hi.3 
t  -  t. 


1+2  -  i 


t  -  t.   N.  .(t) 
i    1,1 

ci+r  ci 


*ftfclMHHfefc*JHIffMdH^ 
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APPENDIX  B 
Subroutine  Index  and  BSPLINE  Program  Call tree 


-83- 


SUBROUTINE  INDEX  FOR  THE  PROGRAM  BSPLINE 

Note:  BOLD-FACE  subroutines  were  written  by  the  author. 

ACTMDB:  subroutine  used  to  turn  on  the  active  mesh  in  MOVING  POINTER  in 
the  MODIFY  BSPLINE  module  of  the  program. 

AIJVAL:  subroutine  to  evaluate  the  reduced  set  of  vertices  needed  in  the 
evaluation  of  derivatives  and  partials  of  the  B-spline  basis 
functions . 

BSPDFM:  subroutine  that  acts  as  the  entry  point  and  control  structure 
for  the  program  module  DISPLAY  FILE  MANAGER. 

BSPGNB:  subroutine  that  acts  as  the  entry  point  and  control  structure 
for  the  program  module  GENERATE  BSPLINE. 

BSPLINE:  main  program  used  to  access  the  other  modules 

BSPMBD:  subroutine  that  acts  as  the  entry  point  and  control  structure 
for  the  program  module  MODIFY  BSPLINE. 

CHGPNT:  subroutine  used  to  graphically  change  the  location  of  control 
points  in  the  command  MOVE  PIJ"S  in  the  module  MODIFY  BSPLINE. 

CMCLOS:  PLOT  10  subroutine  to  end  graphic  communication  and  empty 
graphics  buffer. 

CMDSTR:  subroutine  used  to  parse  the  command  line  for  user  input. 

CMOPEN:  PLOT  10  subroutine  to  open  graphics  communication  with  terminal. 

CMSCMD:  subroutine  to  execute  system  commands  from  FORTRAN  programs. 

CRTDFM:  subroutine  used  to  execute  the  CREATE  FILE  command  in  the 

DISPLAY  FILE  MANAGER  module,  creates  display  arrays  from  the 
existing  program  arrays  if  they  exist. 

DRAW:   PLOT  10  subroutine  to  draw  a  from  present  location  to  location 
specified . 

DRWMSH:  subroutine  to  execute  DRAW  MESH  commands  in  draw  entire  control 
mesh  in  GENERATE  BSPLINE  and  DISPLAY  FILE  MANAGER  modules  of  the 
program. 

DRWPTR:  subroutine  to  draw  the  surface  pointer  in  the  MODIFY  BSPLINE 
module. 


ERSDFM:  subroutine  used  to  execute  command  ERASE  in  the  module  DISPLAY 
FILE  MANAGER  erases  the  screen. 

FRAME:   subroutine  to  draw  the  standard  presentation  window. 
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FRAMER:  subroutine  to  draw  the  full  screen  presentation  window  in  the 
module  DISPLAY  FILE  MANAGER. 

GETDAT:  subroutine  used  to  get  the  counters  used  in  calculating  the 
points  on  the  streamlines. 

GETFIL:  subroutine  used  to  get  filename  and  filetype  when  inputting  and 
outputting  files  to  and  from  the  program 

GETPIJ:  subroutine  that  returns  the  starting  and  ending  control  points 
needed  to  generate  the  given  parameter. 

GETPNT:  subroutine  used  by  the  command  MOVE  PIJ"S  in  the  module  MODIFY 
BSPLINE  to  get  control  point  the  user  wants  to  modify. 

GETSNS:  subroutine  used  to  get  the  sensitivity  used  in  moving  the  points 
in  MOVE  PIJUS  command  in  the  module  MODIFY  BSPLINE. 

GRSTOP:  PLOT  10  subroutine  to  terminate  the  graphics. 

INTDFM:  subroutine  used  to  execute  the  command  RESET  in  the  module 
DISPLAY  FILE  MANAGER,  reintializes  the  display  arrays. 

INTSPL:  subroutine  used  to  execute  the  command  RESET  in  the  module 

GENERATE  BSPLINE,  reinitializes  the  database  and  streamlines  and 
normal  arrays. 

ISOSPL:  subroutine  used  to  execute  the  command  ISOMETRIC  in  the  module 

GENERATE  BSPLINE,  generates  an  isometric  of  the  present  data  set 
by  rotating  first  45  degrees  about  the  y  axis  and  then  35.246 
degrees  about  the  x  axis. 

KAS2AM:  PLOT  10  subroutine  used  to  send  commands  to  the  terminal  from 
the  host. 

LBLMSH:  subroutine  used  by  MOVE  PIJ"S  in  the  module  MODIFY  BSPLINE,  to 
list  the  labels  of  the  active  mesh  in  the  margin  of  the  drawing 
area. 

MDZOOM:  subroutine  used  to  do  a  soft  clip  of  the  drawing  area  or  window. 

MENLST:  subroutine  used  to  draw  the  menu  in  the  presentation  window. 

MNMXER:  subroutine  used  to  find  the  3  dimensional  minimums  and  maximums 
of  the  specified  control  mesh  which  are  then  used  to  size  the 
drawing  window. 

MOVE:    PLOT  10  subroutine  to  move  the  graphics  cursor  without  drawing. 

MPJMDB:  subroutine  used  to  implement  the  command  MOVE  PIJ"S  in  the 
module  MODIFY  BSPLINE. 
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MVPMDB:  subroutine  to  implement  the  MOVING  POINTER  command  in  the  module 
MODIFY  BSPLINE. 

NEWPAG:  PLOT  10  subroutine  to  clear  the  screen. 

NIKVAL:  subroutine  to  evaluate  the  B-spline  basis  function. 

NMLGNB:  subroutine  to  execute  the  command  DRAW  NORMALS  in  the  module 
GENERATE  BSPLINE,  draws  the  normals  to  the  surface. 

NRMVCT:  subroutine  to  calculate  the  normal  vector  given  the  point  on  the 
surface  and  the  parametric  partial  derivatives. 

NUMBER:  subroutine  to  convert  a  character  string  from  the  command  line 
into  a  real  or  integer  number  with  error  checking. 

PRTFIL:  subroutine  to  execute  the  command  PRINT  FILE  which  will  print 
all  the  program  variables  to  the  screen  or  to  a  file. 

RDRDFM:  subroutine  to  execute  the  command  READ  FILE  in  the  module 

DISPLAY  FILE  MANAGER,  inputs  a  display  file  into  the  display 
arrays . 

RDRSPL:  subroutine  to  execute  the  command  READ  FILE  in  the  module 
GENERATE  BSPLINE,  inputs  a  data  file  to  the  program. 

ROTPTS:  subroutine  to  rotate  an  array  of  points  using  a  homogenous 
transformation. 

ROTSPL:  subroutine  to  execute  the  command  ROTATE  in  the  module  GENERATE 
BSPLINE,  rotates  the  data  file  by  a  specified  angle  and  about  a 
specified  axis. 

SAVDFM:  subroutine  to  execute  the  command  SAVE  FILE  in  the  module 

DISPLAY  FILE  MANAGER,  outputs  a  display  file  to  a  specified  file 
on  the  disk. 

SAVSPL:  subroutine  to  execute  SAVE  command  in  the  module  GENERATE 

BSPLINE,  outputs  a  data  file  to  a  specified  file  on  the  disk. 

SCLDFM:  subroutine  to  execute  SCALE  OUTPUT  command  in  the  module  DISPLAY 
FILE  MANAGER,  scales  the  size  of  the  full  screen  output  window 
by  a  factor  between  0.01  and  1. 

SCLPTS:  subroutine  to  scale  an  array  of  points  using  a  homogenous 
scaling  transformation. 

SCLSPL:  subroutine  to  execute  the  command  SCALE  DATA  in  the  module 

GENERATE  BSPLINE,  scales  the  data  file  by  a  specified  amount. 

SCHDFM:  subroutine  to  execute  the  command   SCREEN  TYPE  in  the  module 
DISPLAY  FILE  MANAGER,  used  to  specify  whether  the  standard 
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presentation  window  or  full  screen  presentation  window  is  used 
in  the  module. 

SETSCH:  subroutine  used  to  set  up  the  window  and  viewport  by  finding  the 
minimums  and  maximums  of  the  control  mesh  and  specifying  a  hard 
clipped  window. 

SETUPS:  subroutine  to  put  the  cursor  key  functions  in  the  presentation 
window  in  the  MODIFY  BSPLINE  module. 

SPLPRT:  subroutine  to  evaluate  the  parametric  partial  derivatives  used 
to  calculate  the  surface  normals. 

SPLSFC:  subroutine  to  evaluate  the  point  on  a  surface. 

STBDFM:  subroutine  to  execute  the  command  DRAW  STREAMLINES  in  the  module 
DISPLAY  FILE  MANAGER,  draws  the  surface  streamlines  in  the 
display  file. 

STRGNB:  subroutine  to  execute  the  command  DRAW  STREAMLINES  in  the  module 
GENERATE  BSPLINE,  draws  the  streamlines  generated  using  the 
STREAMLINES  AND  NORMALS  command. 

STKMDB:  subroutine  to  execute  the  command  DRAW  STREAMLINES  in  the  module 
MODIFY  BSPLINE,  draw  the  streamlines  generated  when  modifying 
the  surface. 

STRNML:  subroutine  to  execute  the  command  STREAMLINES  AND  NORMALS  in  the 
module  GENERATE  BSPLINE,  evaluates  the  points  and  normals  to  the 
surface  used  in  display. 

TEXT:   PLOT  10  subroutine  to  put  text  on  the  graphics  screen. 

TRMCHR:  subroutine  to  execute  commands  to  the  terminal  with  the  command 
specified  by  a  character  string. 

TRMCMD:  subroutine  to  execute  commands  to  the  terminal  with  the  command 
specified  by  an  integer  number,  replaced  by  TRMCHR. 

TRML1:   subroutine  to  change  the  terminal  from  4014  graphics  mode  to 
ANSI  mode. 

TRML2:   subroutine  to  change  the  terminal  from  ANSI  mode  to  4014 
graphics  mode. 

TRML3:   subroutine  to  clear  the  entire  ANSI  screen. 

TRML4:   subroutine  to  position  the  ANSI  cursor  at  a  specified  location 
on  the  screen. 

TRML5:   subroutine  to  clear  the  graphics  screen. 
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TRML7:   subroutine  to  erase  the  ANSI  screen  above  the  current  cursor 
location. 

TRML8:   subroutine  to  change  the  graphics  drawing  cursor  to  the  "CLEAR" 
submode. 

TRML9:   subroutine  to  change  the  graphics  drawing  cursor  to  the  "XOR" 
submode. 

TRML10:  subroutine  to  change  the  graphics  drawing  cursor  to  the  "OR" 
submode. 

TRKLll:  subroutine  to  set  the  ANSI  scrolling  region  to  a  specified  area. 

TRML12:  subroutine  to  erase  the  ANSI  screen  below  the  current  cursor 

location. 
TKML13:  subroutine  to  set  the  ANSI  relative  origin  mode. 

TRML14:  subroutine  used  to  set  the  ANSI  absolute  origin  mode. 

TRML15:  subroutine  used  to  set  the  graphics  line  type  in  the  standard 
submode. 

TRML16:  subroutine  used  to  set  the  graphics  line  type  in  the  write 
through  submode. 

TXICUR:  PLOT  10  subroutine  used  to  position  graphics  text. 

TXSIZE:  PLOT  10  subroutine  used  to  specify  the  size  of  graphics  text. 

UPDMDB:  subroutine  used  to  execute  the  UPDATE  command  in  the  module 

MODIFY  BSPLINE,  recalculates  the  streamlines  and  normals  after 
changes  have  been  made  in  MOVE  PIJ"S. 

VWPORT:  PLOT  10  subroutine  used  to  specify  the  portion  of  the  screen  to 
be  used  for  the  window. 

WINDOW:  PLOT  10  subroutine  used  to  specify  the  dimension  of  the 
viewport. 
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CALLTREE  FOR  BSPLINE  PROGRAM 

UTILITIES 

System  Utilities 

CMSCMD ( CMDCMS , I RC ) 

PLOT  10  Subroutines 

GRSTRT( IDEVIC , IOPT) 

GRSTOP 

CMOPEN 

CMCLOS 

'NEWPAG 

VWPORT(XMIN,XMAX,YMIN,YMAX) 

WINDOW(XMINsXMAX  ^MINjYMAX) 

MOVE(PX.PY) 

DRAW(PX,PY) 

TEXTClLENSTjTCFRAY) 

TXSIZE(ISIZE,PXSIZE,PYSIZE) 

TXICUR(IPOS) 

KAS2AM(  ICNT ,  IASRAY  ,  IAMRAY) 

Terminal  Utilities 

TKMCHR(CHRCMD) 
TRMCMD(NCMND) 
TRML1 

Plot 10:    CMOPEN ,KAS2 AM 
TRML2 

Plot 10:    CMOPEN .KAS2AM 
TRML3 

PlotlO:    CMOPEN ,KAS2AM 
TRML4( LNE1 , LNE2 , LNE3 , LNE3 ) 

PlotlO:    CMOPEN ,KAS2AM 
TRML5 

PlotlO:    CMOPEN ,KAS2AM 
TRML7 

PlotlO:    CMOPEN ,KAS2AM 
TRML8 

PlotlO:    CMOPEN ,KAS2 AM 
TRML9 

PlotlO:    CMOPEN ,KAS2AM 
TRMLIO 

PlotlO:    CMOPEN, KAS2AM 
TRML1KTOP1  ,TOP2,BOTl  ,BOT2) 

PlotlO:    CMOPEN .KAS2AM 
TRML13 

PlotlO:    CMOPEN ,KAS2AM 
TRML14 

PlotlO:    CMOPEN ,KAS2AM 
TRML15C NUMBER) 
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PlotlO:    CMOPEN,KAS2AM 
TKML1 6 (NUMBER) 

PlotlO:    CMOPEN.KAS2AM 

Screen  Utilities 

SETSCN(PTS,  NIP,  NJP,  IDMPTS,  SCRN,VWSCL) 

Plot    10:    CMOPEN,CMCLOS,VWPORT, WINDOW 
MNMXER(PTS , I START , I STOP , J START , J STOP , IDMPTS , 

XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX) 
MDZOOM(XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX, SCREEN, SCL, BORDER) 

Plot   10:    CMOPEN,CMCLOS,VWPORT, WINDOW 
MENLST 

Plot    10:    CMOPEN,CMCLOS,TXSIZE, TEXT, MOVE, DRAW, WINDOW, VWPORT 
FRAME 

Plot    10:    CMOPEN,CMCLOS,NEWP AG, MOVE, DRAW, WINDOW, VWPORT 

Terminal:    TBMCMD 
FRAMER(SCN,SCL) 

Plot    10:    CMOPEN,CMCLOS,NEWPAG,MOVE, DRAW, WINDOW, VWPORT 

Terminal:  TBMCMD 

Command  Line  Utilities 

CMDSTR( IPCMND , LCM , CMSTR , ERR) 
NUMBER( NTYPE , NSTR , I , R, NERR) 
GETFIL(FNAME ,FTYPE ,TRYAGN) 
Command:  CMDSTR 

Matrix  Utilities 

ROTPTS(PTS , NIP ,NJP , IDMPTS , ANGLE , AXIS ,  RPTS) 
SCLPTS(PTS, NIP, NJP, IDMPTS, SCLFCT,SPTS) 
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BSPLINE  MAIN 

Plot    10:    GRSTRT,GRSTOP,CMOPEN,CMCLOS 
Terminal:    TRMCHR 
Command:    CMDSTR 
Screen:   FRAME 

Others:    INTSPL,INTDFM  will  be   specified    in  GENERATE  BSPLINE  and 
DISPLAY   FILE  MANAGER 

<GNB>  GENERATE   BSPLINE:    BSPGNB 

<DFM>  DISPLAY  FILE  MANAGER:    BSPDFM 

<PF>  PRINT  FILE:    PRTFIL 
System:    CMSCMD 

Plot   10:    CMOPEN,CMCLOS,NEWPAG 
Command:    CMDSTR 
Screen:    FRAME  ,FRAMER 

<HP>  HELP:    HLPCMD      NOTE:    command   not    implemented 

<QT>  QUIT  PROGRAM 
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GENERATE  BSPLINE:  BSPGNB 

Plot  10:  CM0PEN,CMCL0S 
Terminal:  TRMCHR 
Command :  CMDSTR 
Screen:  SETSCN.FRAME 

<RF>  READ  FILE  RDRSPL 
System:  CMSCMD 
Plot  10:  CM0PEN,CMCL0S 
Terminal:  TRMCMD 
Command:  CMDSTR, GETFIL 
Others:  INTSPL 

<DM>  DRAW  MESH:  DRWMSH(PTS ,NIP ,NJP ,IDMPTS) 
Plot  10:  CMOPEN.CMCLOS, MOVE, DRAW 
Terminal:  TRMCMD 

<DN>  DRAW  NORMALS:  NMLGNB 

Plot  10:  CM0PEN,CMCL0S, MOVE, DRAW 
Terminal:  TRMCMD 

<DS>  :>IAW  STREAMLINES:  STRGNB 

Plot  10:  CMOPEN,CMCLOS, MOVE, DRAW 
Terminal:  TRMCMD 

<IS>  ISOMETRIC:  ISOSPL 

Plot  10:  CM0PEN,CMCL0S 

Terminal:  TRMCMD 

Screen:  FRAME, SETSCN 

Matrix:  ROTPTS 

Others:  DRWMSH(PTS,NIP,NJP,IDMPTS) 

<RT>  ROTATE:  ROTSPL 

Plot  10:  CMOPEN.CMCLOS 

Terminal:  TRMCMD 

Command:  CMDSTR 

Screen:  FRAME, SETSCN 

Matrix:  ROTPTS 

Others:  DRWMSH(PTS,NIP,NJP ,IDMPTS) 

<SD>  SCALE  DATA:  SCLPTS 

Plot  10:  CMOPEN,CMCLOS 

Terminal:  TRMCMD 

Command:  CMDSTR, NUMBER 

Screen:  FRAME, SETSCN 

Matrix:  SCLPTS 

Others:  DRWMSH(PTS,NIP,NJP,IDMPTS) 

<SN>  STREAMLINES  AND  NORMALS:  STRNML 
Plot  10:  CM0PEN,CMCL0S 
Others: 

GETDAT(KNTIMX ,KNTJMX , INTVAL ,MXSTRM ,MXPNTS , 
I STRM , IP NTS , J STRM , JP NTS ) 
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Plot    10:    CM0PEN,CMCL0S 
Terminal:    TBMCMD 
Command:    CMDSTR, NUMBER 

SPLSFCCUI.VJ.PIJ.IORDER.KNOTI.INUMjJORDERjKNOTJ.JNUM.MXNUM.IDIM.S) 

NTKVAL(  T  ,KORDER ,  IKNOT  ,KNOT  ,KNTMAX  ,  JDRV ,  NIK ) 
SPLPRT(UI,VJ,PIJ,IORJ,IORDER,KNOTI,INUM, 

JORDER  ,KNOT J  ,  JNUM  ,MXNUM ,  NDRV ,  ID IM ,  SP ) 
AIJVAL(KORDER,  IKNOT  ,KNOT ,  AI JVRT ,  IDIM ,  JDRV ,    I  J) 
NIKVAL( T  ,K ORDER ,  IKNOT  ,KNOT  ,KNTMAX  ,  JDRV , NIK) 
NRMVCT( UVP NT , PU , P V , ID IM , VCTNML) 

<MDB>  MODIFY   B SPLINE:    BSPMDB 

<PF>  PRINT  FILE:    PRTFIL  same   as    in  BSPLINE  MAIN 

<ER>  ERASE 

Screen:  FRAME, SETSCN 

<RS>  RESET:  INTSPL 

Plot  10:  CMOPEN,CMCLOS 
Terminal:  TRMCMD 

<SV>  SAVE  FILE:  SAVSPL 
System:  CMSCMD 
Plot  10:  CMOPEN,CMCLOS 
Terminal:  TRMCMD 
Command:  CMDSTR, GETFIL 

<QT>  QUIT 

Terminal:  TRMCHR 
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MODIFY  BSLPLINE:  BSPMDB 
Terminal:  TRMCHR 
Command:  CMDSTR 
Screen:  FRAME, SETSCN 
Others: 
SETOPS(NUMSET) 

Plot  10:  CMOPEN.CMCLOS.TEXTjTXSIZE.WINDOW.VWPORT 

Terminal:  TRMCHR 

<MP>  MOVING  POINTER:  MVPMDB( ISTART, ISTOP ,DIVI , JSTART, JSTOP ,DIVJ) 
Terminal:  TRMCHR, TRMCMD 
Command:  CMDSTR 
Moving  Pointer  Utilities 

GETPIJ(T,KNOT,KORDER,NVERT,BGNPIJ,ENDPIJ) 
DRWPTR( UVX , UVY , U VZ , NVX , NVY , NVZ ) 

Plot  10:  CMOPEN,CMCLOS, MOVE, DRAW 

Terminal:  TRMCHR 
Others: 
MSHMDB( PTS, I START, I STOP, JSTART, JSTOP, IDMPTS) 

Plot  10:  CMOPEN,CMCLOS, MOVE, DRAW 

Terminal:  TRMCHR 

<DA>  DRAW  ACTIVE  MESH:  ACTMBD 
Terminal:  TRMCHR 
Command:  CMDSTR 

<DM>  DRAW  MESH:  MSHMDB( PTS , I START, I STOP , JSTART, JSTOP , IDMPTS) 
Plot  10:  CMOPEN.CMCLOS, MOVE, DRAW 
Terminal:  TRMCHR 

<DS>  DRAW  STREAMLINES:  STRMDB(ISTRTS,ISTOPS, JSTRTS, JSTOPS, 

ISTRTP ,ISTOPP , JSTRTP , JSTOPP , IDIM) 
Plot  10:  CMOPEN,CMCLOS, MOVE, DRAW 
Terminal:  TRMCMD 

<MV>  MOVE  PIJ"S:  MPJMDB 

Plot  10:  CMOPEN,CMCLOS 

Terminal:  TRMCHR 

Command:  CMDSTR 

Screen:    FRAME ,MNMXER,MD ZOOM 

Move  PIJ"S  Utilities 

GETSNS(SNS) 

Terminal:  TRMCHR 

Command:  CMDSTR, NUMBER 
GETPNT( ISTART , I STOP , JSTART , JSTOP , ILOC , JLOC ) 

Terminal:  TRMCHR 

Command:    CMDSTR, NUMBER 
CHGPNT(PNTS,IPNT,JP NT, IDIM, ISTART, ISTOP, JSTART, JSTOP) 

Plot    10:    CMOPEN,CMCLOS, MOVE, DRAW 

Terminal:    TRMCHR 
LBLMSH( IBGN , IEND , JBGN , JEND ) 

Plot    10:    CMOPEN,CMCLOS, MOVE, DRAW, TEXT, TXSIZE ,TXICUR 

Terminal:    TRMCHR 
DRWPTR( UVX  , UVY , UVZ , NVX , NVY , NVZ ) 

Plot    10:    CMOPEN.CMCLOS, MOVE, DRAW 
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Terminal:    TRMCHR 
GETPIJ(T,KNOT,KORDER,NVERT,BGNPIJ,ENDPIJ) 
MSHMBD(PTS , I START , I STOP , JSTART , JSTOP , IDMPTS) 

Plot    10:    CM0PEN,CMCL0S, MOVE, DRAW 

Terminal:    TRMCHR 
SETUPS(NUMSET) 

Plot    10:    CMOPEN.CMCLOS, TEXT, TXSIZE, WINDOW ,VWPORT 

Terminal:  TRMCHR 

<UP>  UPDATE:    UPDMBD(lSTRTS,ISTOPS, JSTRTS, JSTOPS, 

ISTRTP ,ISTOPP , JSTRTP , JSTOPP , 
IDVSTR , IDVPTS , JDVSTR , JDVPTS ) 
Plot   10:    CMOPEN.CMCLOS 
Terminal:    TRMCHR 
Others: 
SPLSFCCUI.VJ.PIJ.IORDER^NOTI.INUMjJORDER.KNOTJ.JNUM.MXNUMjIDIM.S) 

NIKVAL(T,KORDER,  IKNOT  ,KNOT,KNTMAX,  JDRV, NIK) 
SPLPRTC DI , VJ  ,PIJ ,  IORJ ,  IORDER ,KNOTI ,  INUM , 

JORDER,KNOTJ ,  JNUM  ,MXNUM ,  NDRV ,  IDIM ,  SP) 
AIJVAL(K  ORDER,  IKNOT  ,KNOT ,  AIJVRT ,  IDIM ,  JDRV ,  AIJ) 
NIKVALC  T  ,K  ORDER ,  IKNOT  ,KNOT  ,KNTMAX ,  JDRV ,  NIK  ) 
NRMVCT( UVP NT , PU , P V , ID IM , VCTNML ) 

<PF>  PRINT  FILE:    same   as    in  BSPLINE  MAIN 

<ER>  ERASE 

Terminal:  TRMCHR 
Screen:  FRAME, SETSCN 

<QT>  QUIT 

Terminal:  TRMCHR 
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DISPLAY  FILE  MANAGER:  BSPDFM 
Plot  10:  CM0PEN,CMCL0S 
Terminal:  TRMCHR, TRMCMD 
Command:  CMDSTR 
Screen:  FRAME, SETSCN 

<CR>  CREATE  FILE:  CRTDFM 

Plot  10:  CMOPEN.CMCLOS 
Terminal:  TRMCHR 

<RF>  READ  FILE:  KDRDFN 
System:  CMSCMD 
Plot  10:  CMOPEN.CMCLOS 
Terminal:  TRMCMD 
Command:  CMDSTR, GETFIL 
Other: 
INTDFM 

Plot  10:  CMOPEN.CMCLOS 

<DM>  DRAW  MESH:  DRWMSH(PTS ,NIP ,NJP ,IDMPTS) 
Plot  10:  CM0PEN,CMCL0S, MOVE, DRAW 
Terminal:  TRMCMD 

<DS>  DRAW  STREAMLINES  STRDFM 

Plot  10:  CMOPEN.CMCLOS, MOVE, DRAW 
Terminal:  TRMCMD 

<ST>  SCREEN  TYPE:  SCNDFM 

Plot  10:  CM0PEN,CMCL0S 
Terminal:  TRMCMD 
Command:  CMDSTR 
Other: 
ERSDFM: 

Plot  10:  CM0PEN,CMCL0S 

Terminal:  TRMCMD 

Screen:  FRAME, FRAMER, SETSCN 

<S0>  SCALE  OUTPUT:  SCLDFM 

Plot  10:  CM0PEN,CMCL0S 
Terminal:  TRMCMD 
Command:  CMDSTR, NUMBER 
Other: 
ERSDFM: 

Plot  10:  CM0PEN,CMCL0S 

Terminal:  TRMCMD 

Screen:  FRAME, FRAMER, SETSCN 

<PF>  PRINT  FILE:  same  as  in  BSLPLINE  MAIN 

<ER>  ERASE:  ERSDFM 

Plot  10:  CMOPEN.CMCLOS 

Terminal:    TRMCMD 

Screen:    FRAME, FRAMER, SETSCN 

<RS>  RESET:    INTDFM 
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Plot  10:  CM0PEN,CMCL0S 
Terminal :  TBMCMD 

<SV>  SAVE:  SAVDFM 

System:  CMSCMD 
Plot  10:  CM0PEN,CMCL0S 
Terminal:  TKMCMD 
Command:  CMDSTR,GETFIL 

<QT>  QUIT: 

Terminal:  TRMCHR 
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APPENDIX  C 
Bell  Data  File 

a.  PRINT  FILE  Output 

b.  Data  file 

c.  Plot  of  I  Parameter 

d.  Plot  of  J  Parameter 
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APPENDIX  C:  display  of  bell  data  file 
for  parameter  J 
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ABSTRACT 

The  objective  of  this  research  was  to  demonstrate  how  the  B-spline 
approximation  could  be  used  to  generate  and  display  3  dimensional 
surfaces  and  to  develop  a  method  to  locally  modify  the  surface 
representation. 

The  work  resulted  in  an  interactive  program  to  generate,  display,  and 
modify,  integral  or  rational,  non-periodic  B-spline  surfaces  through 
order  10.  The  modification  module  utilized  the  inherent  2  dimensional 
characteristics  of  the  parametric  surface  representation  to  move  a 
pointer  on  the  3  dimensional  surface.  The  location  of  the  pointer  was 
used  to  specify  a  finite  set  of  control  points  which  could  be  used  to 
modify  the  surface  by  changing  their  positions  in  space. 

The  program  ran  on  the  Kansas  State  University  Computing  Center's 
mainframe  computer  using  a  sub  set  of  PLOT  10  IGL  subroutines  and  a 
Selenar  Hirez  100  graphics  terminal  as  the  user  interface. 


